Fire Containment in Grids of Dimension Three and 

Higher 



Mike Develin Stephen G. Hartke* 

American Institute of Mathematics Department of Mathematics 

360 Portage Avenue University of Illinois 

Palo Alto, CA 94306 Urbana, IL 61801 

develin@post . harvard . edu hartkeOmath . uiuc . edu 

September 9, 2004 



Abstract 

We consider a deterministic discrete-time model of fire spread introduced by Hart- 
nell [1995] and the problem of minimizing the number of burnt vertices when deploying 
a limited number of firefighters per timestep. While only two firefighters per timestep 
are needed in the two dimensional lattice to contain any outbreak, we prove a conjec- 
ture of Wang and Moeller [2002] that 2d — 1 firefighters per timestep are needed to 
contain a fire outbreak starting at a single vertex in the d-dimensional square lattice for 
d > 3; we also prove that in the ci-dimensional lattice, d > 3, for each positive integer 
/ there is some outbreak of fire such that / firefighters per timestep are insufficient 
to contain the outbreak. We prove another conjecture of Wang and Moeller that the 
proportion of elements in the three-dimensional grid P n x P n x P n which can be saved 
with one firefighter per time step when an outbreak starts at one vertex goes to as n 
gets large. Finally, we use integer programming to prove results about the minimum 
number of timesteps needed and minimum number of burnt vertices when containing 
a fire outbreak in the two dimensional square lattice with two firefighters per timestep. 

Mathematics Subject Classification: 05C75, 90C35 

Keywords: firefighter, containment strategy, vaccination strategy 



1 Introduction 

Hartnell [5] introduced a deterministic discrete-time model of fire spread on a graph G and 
considered how firefighters can act to stop a fire outbreak. In this model, an outbreak of 
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fire starts at a set of root vertices of G at time t = 0. In response, firefighters are placed 
at the vertices a^i, ai^, • • • , Qi, Cl at time t = 1, where the firefighters defend or protect each 
vertex from the spreading fire. The fire then spreads from burning vertices to non-defended 
neighbors. Firefighters are again deployed to defend the vertices 02,1, . . . , a2, C2 at time t — 2 
(the vertices a^i, 01,2, . . . , ai, Cl remain defended), and the fire spreads again. The process 
continues until the fire can no longer spread. We say that the fire outbreak is contained 
after t time steps if there is some finite time t such that after the disease spreads during 
time t, only a finite number of vertices are burnt and the disease can no longer spread. The 
motivating question is to find an optimal sequence of defended vertices that minimizes the 
total number of burnt vertices. 

The fire spread model is also relevant in epidemiology. Traditionally, epidemiological 
models assume that the population being studied is well-mixed in the sense that any pair of 
individuals are just as likely to come in contact and transmit a disease as any other. However, 
recently epidemiologists have attempted to incorporate spatial information into their models 
[712112]. The model of fire spread presented above can be considered as modeling a perfectly 
contagious disease with no cure, where vertices adjacent to infected vertices become infected 
at every discrete time step and, once infected, remain infected from then on. The response 
allowed is only a limited number of vaccinations of non-infected vertices per time step. The 
limited number of vaccinations is particularly relevant to real-world situations because of 
limited availability of the vaccine or limited numbers of health personnel to administer the 
vaccine. The main question is of course to minimize the total number of infected vertices. 

In this work we study fire containment on square grids. Grids are a natural class of graphs 
on which to consider both fire and disease spread since they are often used to represent 
geographic areas. Both Wang and Moeller [Hj and MacGillivray and Wang [B] studied grids 
to find algorithms for containment. Wang and Moeller showed that two firefighters per time 
step is sufficient to contain a fire outbreak in a two dimensional square grid, and conjectured 
that 2d — 1 firefighters are necessary to contain a fire outbreak in a d dimensional square 
grid. We prove this conjecture in section |21 Fogarty [3] showed that two firefighters suffice in 
the two dimensional square lattice to contain any finite outbreak of fire where an arbitrarily 
large but finite number of vertices are initially on fire. However, we prove that for any fixed 
number / of firefighters, there is a finite outbreak of fire in which / firefighters per time step 
are insufficient to contain the outbreak. We also prove the conjecture of Wang and Moeller 
that the proportion of elements in the three-dimensional grid P n x P n x P n which can be 
saved by using one firefighter per time step when an outbreak at one vertex occurs goes to 
as n gets large. 

In section |3] we provide an alternate proof using integer programming of Wang and 
Moeller's result that the minimum number of time steps needed to contain a fire outbreak 
in a two dimensional square grid when using two firefighters per time step is 8. We also use 
this technique to prove that the minimum number of burned vertices in such an outbreak is 
18. 

We use the following terminology to describe the fire spread and firefighter response. 
During the t th time step for t > 0, firefighters are deployed and then the fire spreads. If 
we describe the state of vertices at the beginning of the t th time step, we mean before the 
firefighters are deployed during the t th time step. If we describe the state of vertices at 
the end of the t th time step, or equivalently, at the end of t time steps, we mean after the 
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fire has spread during the t th time step. A firefighter may defend neither a burnt vertex 
nor a previously defended vertex. Once fire has spread to a vertex v, we say that v is a 
burnt vertex. After being burnt or defended, a vertex remains in that state until the process 
ends. In addition to the burnt and defended vertices, we say that a vertex v is saved at the 
end of the t th time step if there is no path from v to the root consisting only of burnt and 
non-defended vertices at the end of the t th time step. 

We consider the infinite d- dimensional square grids L d . The vertices of h d are the points 
of M. d with integer coordinates, and x is adjacent to y if and only if x is distance 1 from y in 
the usual Euclidean i 2 metric. 

2 Three and Higher Dimensional Square Grids 

Wang and Moeller proved in [H] that an outbreak starting at a single point in a regular 
graph of degree r can be contained with if r — 1 firefighters can be deployed per time step. 
Specifically, for the d dimensional square grid h d , 2d — 1 firefighters suffice to contain an 
outbreak starting at a single point. They conjectured that this bound is tight, and we present 
a proof of this conjecture here. 

Wang and Moeller observed that at least two firefighters per time step are needed for 
containment in L 3 , and Fogarty showed in [3] that at least three firefighters per time step 
are needed to contain the outbreak. Her main theorem involves a "Hall-type condition" for 
the graph, which provides a lower bound for how fast the fire can spread. The theorem 
considers the front of the fire, which is the set of burnt vertices farthest from the root. The 
theorem states that if this front grows quickly (i.e., by at least /) regardless of its precise 
shape, then it cannot be contained by deploying / firefighters per time step. Theorem El 
strengthens Fogarty's theorem by considering initial growth of the fire that is faster than / 
so that the fire reaches a "critical mass" and can sustain growth of the front by at least / 
from that point onward. 

First we state some definitions. 

Definition 1. Let D k denote the set of vertices in a rooted graph G that are distance k 
from the root vertex r. Let r k denote the number of firefighters in D k +i, D k+2 , • • • at the end 
of the k th time step. These firefighters can be thought of as "reserve" firefighters since they 
are not adjacent to the fire when deployed. We define r to be 0. Let B k C D k denote the 
number of burned vertices in D k at the end of the k th time step. 

Theorem 2. Let G be a rooted graph, h a positive integer, and do, a±, . . . , ah positive integers 
each at least f such that the following holds: 

1. Every A C D , A^$, satisfies \N(A) H D x \ > \A\ + a . 

2. Fori <k < h, every A C D k where \A\ > 1 + Eto( a * ~ /) satls fi es \ N ( A ) n D k+i\ > 
\A\ + a k . 

3. For k> h, every A C D k such that \A\>\ + E?=o( a i ~ /) satisfies \N(A) n D k+l \ > 
14+/- 
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Suppose that at most f firefighters per time step are deployed. Then 

{1 ifn = 0, 

l + r n + E"=o(^-f) ifl<n<h+l, (1) 
l + r n + Eto(«<-/) ifn>h + l, 

regardless of the sequence of firefighter placements. Specifically, f firefighters per time step 
are insufficient to contain an outbreak that starts at the root vertex. 

Proof. Let p n +i denote the number of firefighters placed in D n+ \ at time n + 1, and let p< n 
denote the number of reserve firefighters placed in D n+ \ during time steps 1, . . . ,n. Note 
that 

r n +i < (r n ~ P<n) + (/ - Pn+i) = r n + / - p n+1 - p< n . (2) 

This follows since r n — p< n is the number of firefighters placed in D n+ 2, D n+ s, . . . for times 
1, . . . , n, and at most / — p n+ i firefighters are available to be placed in D n+2 , D n+3 , ... at 
time n+ 1. Strict inequality occurs if a firefighter is placed in for k < n + 1 at time n + 1. 

We prove by induction on n. For n — 0, \B \ = 1 holds trivially. We assume the 
result holds for n, < n < h, and prove the result for n + 1. By the inductive hypothesis, 

fl if n = 0, 

\B n \ > < i 3 

\l + rn + Er=o(«i-/) ifl<n</i, 
and so by hypotheses [T] and |2l 

|iV(S n )nD n+ i| > \B (4) 



Thus, 



= \N(B n ) n D n+ i| -p n+1 -p<n 
> \B n \ +a n - p n+1 - p< n , by (gj), 

n-1 

<^n Pn+1 P<ni 

by ©, 



i=0 

n-1 



= 1 + {r n + f - p n+1 - p< n ) + ^(tti - f) + (a„ - /) 

i=0 

n 

> 1 + r n+ i 

+ ^(^-/), by©. 



i=0 



This proves for < n < h + 1 . 

We now prove © for n > h + 1 using induction on n. Note that ([T} holds for n = /i + 1 
from above. We thus assume © holds for n > h + 1, and we prove the result for n + 1. By 
inductive hypothesis, 



h 



\B n \>l + r n + J2( a i-f), ( 5 ) 



i=0 



and so by hypothesis 01 © holds for n > h. Thus, 

\B n+1 \ = \N(B n ) nD n+1 \ -p n+ l ~P<n 

> \B n \ + f - p n+1 - p< n , by (gj), 

h 

> 1 + r n + ^(cti -/) + /- Pn+i - P<«, by ©, 

i=0 

= 1 + (r n + / - p n+ i - p< n ) + ^2(di - /) 

i=0 

h 

>i + r n+1 + J2( a i-f)i h y©- D 

This completes the proof of Theorem El 

We now turn our attention to square lattices of dimension three and higher. It will prove 
convenient to partition these lattices into identical subgraphs. 

Definition 3. The orthants of M. d are the 2 d regions defined by the hyperplanes Xi = —1/2 
in IR d , i = 1, . . . , d. Let the orthants in h d be the subsets of vertices that lie in each orthant of 
]R d . Thus, the j th coordinates of all the vectors in a given orthant of M. d are all non-negative 
or are all negative, for j = 1, . . . , d. Let denote the vertices of C h d in the orthant 
whose elements are all non-negative. 

Let v = (vi,v 2 , ■ ■ ■ , Vd) be an element of D k C L d . Let Cj(-u) denote v iy and for a set 
A C D k define A\ = {v e A : Ci(v) = r}. Let iUj denote (v l , v 2 , . . . , v v i+1 , ...,v d ) G D k+l , 
where v\ — v i + 1 if > or = V{ — 1 if < 0. Thus, is in the same orthant as v. 

Lemma 4. Inh d ford> 3, if A C D k where \A\ > 2d-2, then \N(A) n D k+1 \ > \A\+2d-2. 

Proof. Given any nonempty set A C D k C L ' completely contained in one orthant, we will 
show that 

\N{A) n D k+1 \ >\A\ + d-l, for any d. (6) 
We form a set 5 C N(A) PI -Dfc+i in the following way: 

1. For each u6 A, add to B. 

2. For each 2 < j < d, let rj be the value of the j th coordinate of elements of A that is 
greatest in absolute value. For each v G A J r ., add v^j to B. 

Each vector added to B in step 1 is unique, and each vector added to B in step 2 is also unique 
since the j th coordinate was chosen to be maximum. Thus, \N(A) H D k+ \\ > \B\ > \A\+d—l. 

Let A C Dfc C L d . If A is not completely contained in one orthant, then let A be 
partitioned as 

A = Aj U A 2 U • • • U A g , 
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where each A e is in a different orthant O e . By ©, \N(A e ) H £> fe+ i| > \M + d - 1. Note 
also that the corresponding sets Bi in the proof above for Ai do not overlap since they are 
in different orthants. Hence, 



\N(A) n Ami > ^ \ N ( A i) n °< n 



£=1 

>^P,|+d-l] 

£=1 

> AI + 2d- 2. 



Thus, we may assume that A is completely contained in one orthant, and, without loss of 
generality, we assume that all coordinates of elements of A are non-negative. 

We now proceed to prove the lemma by induction on d. Let A C D£ C h d , where 
\A\ > 2d— 2. Suppose that d = 3. Let rij denote the number of nonempty A\, or, equivalently, 
the number of distinct i th coordinates of elements of A. Let %' be a coordinate where is 
maximized. We claim that n*/ > 3. If ny is 1, then A contains only one element, which 
is a contradiction since \A\ > 2d — 2 = 4. If is 2, then each coordinate has only two 
different values it can assume. However, the sum of the coordinates must remain k. It is 
straightforward to verify that the maximum number of elements in A is 3, which contradicts 
the fact that \A\ > 4. Thus, ?v > 3. 

For each r where A*' is nonempty, form a set A*' C D^lJ C L d_1 by eliminating the 
%' coordinate of each element in A\ ; thus, the A*"s are the parts of A contained in the 



slices of h d taken in direction i'. By (jBJ, iV(A*') PI 



> 



Ai' 



d — 2. For each w in 



N(Af] 



n -D^_J +1 , form an element v in A^(A*') fl Df +1 by inserting r as the %' coordinate. 



Notice that these elements are distinct when the % 
maximum r such that A 



coordinates are distinct. Let m be the 



is nonempty, or equivalently, the largest i' coordinate. For each 
v G A l m , we also have V-+y G N(A) fl -D^+i, and these vectors are distinct from any formed 
above because the %' coordinate is larger. Thus, 



\N(A)nD k+1 \ > 



AI 



> \A\ + n v {d 



d 



AI 



A, 



(7) 



Since |A^| > 1, © implies that 

\N(A) PI D k+ i\ > \A\ +3d- 5, (8) 

and when d = 3, 

|iV(A) n D fc+1 | > |A| + 4 = |A| + 2rf - 2. 

Now suppose that d > 3. Again let denote the number of nonempty A*, and let i' be 
a coordinate where rij is maximized. If ?v > 3, then using the same construction as in the 
d = 3 case, we have ©, and since d > 3, \N(A) n > |A| + 2cf — 2. If = 1, then A 
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contains only one element, which is a contradiction since \A\ > 2d — 2 > 4. We are thus left 
with the case ny = 2. Let m be the maximum r such that A % r is nonempty, or equivalently, 
the largest %' coordinate of elements of A, and let r' ^ m be the minimum value of r where 
A l l is nonempty. If \A l m \ > 2, then using the same construction as in the ny > 3 case, we 
have by © 



A! 



|JV(A)n£> fc+ i| > \A\ + n it {d - 2) + 

> |A| + (2d -4) + 2, since 

> |A| + 2d-2. 



A % „ 



>2, 



If |AjJ = 1, then we again use the construction from the > 3 case. However, 



At 



> 



2d — 3, so by induction, 



N(At) n D 



k-r'+l 



> 









2d — 4. Here, the notation T)t 1 means 



the set D z C L d_1 , emphasizing the dimension of L d_1 . For each v in JV(A*',) fl D^., , l5 
form an element u in iV(^4*,) fl -D^ +1 by inserting r' as the %' coordinate. Additionally, for 
the single vector v 6 and 1 < j < d, v^j G AT (A) fl D^+i, and these vectors are distinct 
from those formed above because the %' coordinate is larger. Thus, 



\N{A)HD 



fc+i 



> 

= \ A \ 
> \A\ 



A„i 



2d -4 



d 



3d 
2d 



3, since 



A„i 



L4| + l, 



2, since d > 3. 



□ 



Lemma |U provides the long-term growth of the front, Condition 3 needed for Theorem 121 
The next lemma gives the complementary requirements. 

Lemma 5. In h d for d>3, if A C D 1 where \A\ > 2, then \N(A) n D 2 \ > \A\ + 4d - 6. 

Proof. Let A C C L d where |A| > 2. Every vector i> £ A is of the form (0, . . . , Xi, . . . , 0), 
where rcj = ±1. Each vector v in A has 2{d — 1) neighbors in Z?2 formed by replacing each 
of the zero coordinates in v with ±1, and one neighbor formed by replacing 1 in the i th 
coordinate with 2 or replacing —1 with —2. If v and v' are vectors of A with nonzero entries 
in different coordinates, then v and v' share exactly one neighbor in D 2 . If v and v' have 
nonzero entries in the same coordinate, then v and v' share no neighbors in D 2 . Thus, 



\N(A) fl D 2 \ > \A\ (2(d 



L4I 2d 



> 1^1 + 1^1 

It is straightforward to use calculus to verify that 

lAIIM-M-?' 

2 2 



1)4 


-1)- 


14 


4) 


2 




(2d 


\A\ 
2 



> 4d - 6, 



7 



where d > 3 and 2 < |A| < 2d, and so 

|N(A) nD 2 | > \A\ + 4d-6. □ 

Theorem 6. In L , 2d — 1 firefighters are needed to contain an outbreak of fire starting at 
a single vertex. 

Proof. Since L d is vertex transitive, we may assume that the root vertex where the fire 
outbreak starts is the origin. We use Theorem [21 with / = 2d — 2, h — 1, ao = 2d — 1, and 
a\ = 4d — 6. The one element set Do has 2c? neighbors in D\ so hypothesis ^ of Theorem [21 
holds, Lemma 03 shows hypothesis El of Theorem [21 holds for k — 1, and Lemma 0] shows 
hypothesis H3 holds for k > 1. By Theorem [21 2d — 2 firefighters are insufficient to contain an 
outbreak starting at the origin. □ 

Fogarty also showed in [3] that two firefighters suffice in L 2 to contain any finite outbreak 
of fire where an arbitrarily large but finite number of vertices are initially on fire. However, 
we prove for h d where d > 3 that for any fixed number / of firefighters, there is a finite 
outbreak of fire in which / firefighters per time step are insufficient to contain the outbreak. 

First we establish the following lemma. Essentially, the lemma says that if we have a 
"front" of x elements, then it will grow outwards by at least Vt{^/x) in the next time step. 

Lemma 7. Let f be any positive integer. If A C D£ C L 3 where \A\ > |(/ — 1)(/ — 2), 
then \N(A) H D+ +1 \ > \A\ + /. 

Proof. Let A C D+ C L 3 be a set where |A| > \{f - l)(f - 2). The elements of B := 
{v^i : v G A} are distinct vertices in N(A) D D£ +1 , and the set B has cardinality equal to 
\A\. Therefore, it suffices to show that if \A\ > |(/ — 1)(/ — 2), then there are at least / 
distinct elements of the form v^j which are not elements of B, where v G A and j G {2, 3}. 

Let m be the largest first coordinate of elements of A, and let t be the smallest first 
coordinate of elements of A. Recall that the sets A\, r = t, t + 1, . . . , m, partition A. Let o~ r 
equal \A].\, so that Y^T=t a r = 1^1- Note that a t ,o~ m > 0. 

Suppose some a r is equal to zero, where t < r < m. Then A is partitioned into the sets 
A\ consisting of all elements of A with first coordinate greater than r and A2 consisting of 
all elements of A with first coordinate less than r. Clearly N(Ai) H N(A2) H D^ +1 = 0. 
Define A[ := {f_>2 : v G Ai} and A' 2 := {^1 : v G ^2}, so that A\ and t4' 2 are subsets of 
D^ +1 . Since is simply a translate of Ai by 1 in the first coordinate, N(A[) fl D^ +2 is a 
translate of iV(Ai) Hi9^" +1 by 1 in the first coordinate. Similarly, N(A' 2 ) C\D^ +2 is a translate 
of N(A2) PI by 1 in the second coordinate. Thus, we have that 

\N(A[ U A' 2 ) n D+ +2 \ < \N(A[) n u+al + \n(a' 2 ) n ^1 

= |iV(A 1 )n J D+ +1 | + |iV(A 2 )n J D+ +1 | 
= |iV(A)n J D+ +1 |, 

where the last equality follows since iV(Ai) fl -D^" +1 and N(A 2 ) H -D^ +1 do not intersect. 
However, A[ U A 2 has the same size as A, but the separation between the largest first 
coordinate of elements of A[ U A' 2 and the smallest first coordinate of A[ U A 2 is less than 
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m — t. Therefore, by induction on m — t we reduce to the case where no oy is equal to zero, 
i.e., there is an element of A with first coordinate r for every t < r < m. 

Consider the sets S r = {tUj : v G Aj,j G {2,3}} C N(A) n D+ +1 . Observe that the 
cardinality of S r is at least oy + 1. Clearly all S r are disjoint, since all elements of S r have 
first coordinate r. The elements of St have t as their first coordinate, while all elements of 
B have first coordinates at least t + 1, so no elements of St are in B. Furthermore, for all 
r > t, if an element of S r is in 5, then by considering its first coordinate, the element must 
be in the set {v^i : v G In particular, this set has size ay_i. If o r + 1 > oy_i, then 

there are at least a T + 1 — o r _i elements in S r not in 5. Therefore, the number of elements 
in N(A) fl Dfr +1 that are not in B is bounded below by 

m 

g(a) : = max (0, o r + 1 — oy_i), (9) 

r=t 

with the convention that a t -i — 0. 

Now take any nonzero sequence a t ,a t+ i, . . . ,a m . We claim that if g(a) < f, then 
J2r=t a r — \U' ~~ !)(/ — 2), which would complete the proof of the theorem. Suppose 
we have some sequence a t , Ot+i, • • • , cr m with g{a) < f. First, suppose that there exists some 
r > t where oy > cr r _i. Then adding 1 to cr r _i decreases the r-th term of (0) by 1, possibly 
adds 1 to the (r — l)-st term, and leaves all other terms unchanged; in particular, it does not 
increase the value of g(a) and increases Yl °r- Therefore, we can reduce to the case where a 
is strictly decreasing. 

Next, suppose we have ay < <J r -i — 1 for some t < r < m. Then adding 1 to oy leaves 
all terms of JHJ) unchanged. Similar to before, this operation does not change g(cr), while 
increasing ^ o r . Doing this repeatedly, we reduce to the case where 

oy„i = oy + 1 (10) 

for all t < r < m. However, this case is easy to evaluate; each term in Q is zero except the 
r = t term, which is equal to a t + 1. Since g(a) = a t + 1 < /, a t < f — 1. Since a m > 0, 
YlT=t °r ^ S a ^ mos t the sum of the first / — 2 positive integers. Thus, 

r=t 

This allows us to prove the following theorem. 

Theorem 8. For any dimension d > 3 and any fixed positive integer f , f firefighters per 
time step are not sufficient to contain all finite outbreaks in L d . 

Proof. Since L 3 is contained in h d for d > 3, it suffices to prove the statement for d = 3. 
We consider an initial outbreak consisting of all of D~£ for k large enough so that \D~£\ > 
|(/ — — 2). To show that / firefighters are insufficient to contain this outbreak, we 
will construct a related graph that captures the essential disease dynamics and then invoke 
Theorem |2l Let G be the subgraph of L 3 induced by vertices with non- negative coordinates 
that are distance at least k from the origin. Let G' be the graph formed from G by identifying 
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all of the vertices in as a single vertex r. An edge exists between vertices x and y in G' 
if xy is an edge in G or if x = r and y e Ng(D^). Let _D 4 ' denote the set of vertices in G' 
that are distance z from the root r. By Lemma 

|# W) n > \d+\ +f>\(f- - 2) + /, 

and so 

|JV(r)nl3i|>(|^ |-l) + i(/-l)(/-2) + /. 

If A' C D^, where i > and \A'\ > -(/ - 1)(/ - 2), then A' corresponds to a set A C _D+_. 
and by Lemma [3 

|A^)nD+ +m |>L4| + /, 

and hence 

|iV(A')n^ +1 | > L4'| + /- 

By Theorem |2] with /i = 0, and Oq = \{f — 1)(/ — 2) + /, / firefighters are insufficient to 
contain an outbreak starting at r in G", and hence / firefighters are insufficient to contain 
an outbreak consisting of all of D\ in L 3 . □ 

The essential problem here is that for d > 3, the boundary of an outbreak grows faster 
than the constant number of firefighters deployed at a given time step. Indeed, in dimension 
d, the boundary grows as a polynomial of degree d— 2. This motivates the following ambitious 
conjecture. 

Conjecture 9. Suppose that f(t) is a function on N with the property that goes to as 
t gets large. Then there exists some outbreak on L d which cannot be contained by deploying 
f(t) firefighters at time t. 

A weaker conjecture would require f(t) to be a polynomial. 

Lemma [7| also allows us to resolve another conjecture of Wang and Moeller in jS]. They 
conjectured that as n gets large, the proportion of elements in the three-dimensional grid 
P n x P n x P n which can be saved by using one firefighter per time step when an outbreak at 
one vertex occurs goes to 0. We prove this conjecture in the following 

Theorem 10. Let v be any vertex of P n x P n x P n , for n > 1. Then the maximum number of 
vertices which can be saved by deploying one firefighter per time step with an initial outbreak 
at v grows at most as 0(n 2 ). In particular, the proportion of vertices which can be saved 
goes to as n gets large. 

Proof. We prove the theorem in the case v = (0,0,0). The general statement easily follows 
by splitting P n x P n x P n into orthants with apex v. We actually prove a stronger statement. 
Consider the graph G induced from the lattice L 3 by vertices with non-negative coordinates 
and distance at most 3n from the origin v. We prove the theorem for the graph G. Note 
that G contains P n x P n x P n as an induced subgraph. 

We claim that \B t \ — r t > * + 2 f+2 for all t regardless of what firefighter placements are 
made. Since there are (*^ 2 ) = f2+3f+2 vertices in Df, this statement is saying that at the end 
of the t th time step the number of reserve firefighters together with the unburned vertices 
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(including defended vertices) in Df cannot exceed t. By considering time up to t = 3n, 
when all vertices have had a chance to be burned, at most 1 + 2 + . . . + 3n = 0(n 2 ) vertices 
are unburned. This implies the same conclusion for P n x P n x P n . 

The proof of the claim is by induction. At the end of the th time step, there are no 
reserve firefighters, and one vertex in D\ is burned; the difference isl — = 1 > 1 = + 2 0+2 
as desired. 

Suppose t > 0, and suppose that the statement is true for t. Then 

t 2 4- t + 2 1 

\B t \> t -^l±^>U(t + l). (11) 

Let / = t + 2. By Lemma 

\N(B t )nD+ +1 \ > \B t \ + f. (12) 

As in the proof of Theorem El let p t +i denote the number of firefighters placed in Df +1 at 
time t + 1, and let p< t denote the number of reserve firefighters placed in DJ, X during time 
steps 1, . . . , t. Thus, 

\B t+l \ - r m = [\N(B t ) n D+ +1 \ - Pt+1 - p< t ] - r t+1 

> \N{B t )nD+ +1 \-r t -l, by0, 

> |5 t |+/-r t -l, by (HI, 

+2 _|_ + i 2 

> ~ 2 + (t + 2)-l, by (HU), 



> 



(t+ 1) 2 + (t + 1) + 2 



Hence the claim follows. □ 

In practice, one can ensure when an outbreak starts at (0, 0, 0) that t vertices in Df are 
unburned at time t. However, because the fire doubles back on itself, it is unclear that one 
can actually save a quadratic number of vertices. Wang and Moeller exhibit the construction 
of building a "fire wall" by defending all of the vertices at distance k from (n, n, n). In order 
for this to be effective, we must be able to cover all ( fc+1 H fc + 2 ) suc ] 1 vertices in the 3n — k 
time steps it takes the fire to reach this hyperplane. This yields k = 0{^Jn). The number 
of vertices saved is the number of vertices at distance k or less from (n,n,n), which is 
(fc+i)(fc+2)(fc+3) ^ fjkjg j s (3(^3) _ 0(n 3 / 2 ). Therefore, the optimal number of vertices saved 
given an initial outbreak at (0, 0, 0) in the grid graph P n x P n x P n when deploying one 
firefighter per time step is between 0(n 3 / 2 ) and 0{n 2 ). 



3 Two Dimensional Square Grid 

According to Wang and Moeller in jH], Hartnell, Finbow, and Schmeisser first proved that 
an outbreak of fire in L 2 starting at a single vertex can be contained using two firefighters 
per time step. Their sequence of firefighter placements contained the outbreak at the end 
of 11 time steps. Wang and Moeller showed that the disease cannot be contained at the 
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end of 7 time steps when using two firefighters per time step and presented a sequence of 
firefighter placements that attains this minimum. Their sequence allows 18 vertices to be 
burned. Surprisingly, Wang and Moeller do not comment on whether their solution attains 
the minimum number of burned vertices. In fact, 18 is the minimum number of burned 
vertices, and we prove this using integer programming. The same technique also gives a 
computer proof of Wang and Moeller's result that at least 8 time steps are needed. Their 
proof relies heavily on case analysis. 

The tightness in the following theorem is due to Wang and Moeller [S]. 

Theorem 11. In h 2 , if an outbreak of fire starts an a single vertex, then when using two 
firefighters per time step at least 18 vertices are burned. This bound is tight. 

Proof. We formulate an integer program using the boolean variables b x t and d X)t . The 
variable b Xs t is 1 if and only if vertex x is burned at or before time t, and d Xi t is 1 if and 
only if x is defended at or before time t. We wish to minimize the total number of vertices 
that become burned. For the integer program to be implementable with a finite number of 
variables and constraints, we restrict the graph to L = {(x, y) G L 2 : \x\ < i and \y\ < £} and 
< t < T, where i and T are chosen to be sufficiently large that the fire never reaches the 
boundary and is completely contained by time T. In the actual computations performed, 
t = 6 and T = 9 proved sufficient. We choose T > 8 to ensure that the fire is actually 
contained and does not grow in the last time step. 
The integer program is 

minimize b Xt T 

subject to:b x>t + d x>t — b Vit _i > 0, for all x G L, y G N(x), and 1 < t < T, (13) 

b x ,t + d X:t < 1, for all x G L and 1 < t < T, (14) 

b x ,t — b x j~i > 0, for all x G L and 1 < t < T, (15) 

dx,t — dx,t-i > 0, for all x G L and 1 < t < T, (16) 

(dx,t - 4,t-i) < 2, for 1 < t < T, (17) 



6 Xj0 = { ' ^ f° r a h ^ G L, (18) 

otherwise, 

d x ,o = 0, for all x G L, (19) 

b x ,t, d Xjt G {0, 1}, for all x G L and < t < T. (20) 

Condition (fTB^) enforces the spread of the fire while respecting vertices defended by a fire- 
fighter. Note that vertices can spontaneously combust, catching fire, but the minimization 
of the objective function ensures that this does not happen in the optimal solution. Condi- 
tion (|14j) prevents a firefighter from defending a burnt vertex, while conditions (|15j) and (|16j) 
ensure that once a vertex is burnt or defended, it stays in that state. Condition (|17|) only 
allows two firefighters per time step. Conditions (fTHj) and ()19|) give the initial conditions at 
time t = 0, and condition (J2*U|) makes the program a binary integer program. 
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Figure 1: Optimal solution of the integer program used in the proof of Theorem fTTl The 
fire outbreak starts at time at the root, and then spreads to the black vertices at the times 
written next to the vertices. The square firefighters Oj are placed at time i. This placement 
of two firefighters per time step in L 2 completely contains the outbreak in 8 time steps, 
allowing only the minimum number of 18 burned vertices. 



The integer program was solved in about 1.83 hours using the GNU Linear Programming 
Kit jlj running on a Pentium IV 2.6GHz processor, and 18 was the minimum number of 
burnt vertices at time t — 9. Figure ^ shows the minimum solution. The fire was completely 
contained and thus did not reach the sides of L. Also note that the solution presented by 
Wang and Moeller in jS] also allows only 18 burnt vertices but is slightly different from the 
solution presented here. □ 

Lemma 12. If an outbreak of fire in L 2 is contained by 14 defended vertices and (x,y) is a 
burnt vertex, then \x\ < 5 and \y\ < 5. 

Proof. Suppose that (x,y) is a burnt vertex, and, without loss of generality, that x > 5. 
Since (x, y) is burnt, there is a path v = (x, y),V\, v 2 , ■ ■ ■ , v t = (0, 0) from (x, y) to the origin 
consisting of burnt vertices. For each < a < 6, there is a vertex v p t a \ such that the first 
coordinate of v p ( a ) is a. Since the fire is contained, there must be a defended vertex above 
and below each of these seven vertices, and there must be at least one defended vertex with 
first coordinate less than and one with first coordinate greater than x. But this requires 
16 defended vertices, resulting in a contradiction. □ 

Theorem 13 (Wang and Moeller). In L 2 , if an outbreak of fire starts at a single vertex, 
then the fire cannot be contained at the end of 7 time steps when using two firefighters per 
time step. Thus, at least 8 time steps are needed to contain the fire, and this bound is tight. 

Proof. We use a similar integer program to the one used in the proof of Theorem ^2 By 
Lemma El if the outbreak can be contained after 7 time steps, then no burnt vertex will 
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have either coordinate equaling 6 in absolute value. We thus use the finite grid L where 
£ = 6, and we use the objective function 



minimize 



x=(a,b)£L 
|a|=6 or |6|=6 



If the disease can be contained after 7 time steps, then the optimal value of the objective 
function will be 0. All of the conditions from the previous integer program are included 
except condition (JT7j) is changed to 



This prevents firefighters from being used after 7 time steps. 

The integer program with T = 9 was solved in about 40 minutes using the GNU Linear 
Programming Kit running on a Pentium M 900MHz processor. The minimum value was 1, 
meaning that in every feasible solution, the fire burned a vertex with one coordinate equaling 
6 in absolute value. This contradicts Lemma and so at least 8 time steps are needed to 
contain an outbreak in L 2 when using two firefighters per time step. □ 



There are many avenues for future work in models of responses to fire and disease spread. 
For infinite graphs, we can ask the same question as for the infinite square grids: What is 
the minimum number of firefighters needed per time step so that only a finite number of 
vertices are burned? Percolation is a related topic whose methods may also apply here. 

From the viewpoint of an arsonist or bioterrorist, one would like to find the most vul- 
nerable vertices in a graph G. A vertex v is most vulnerable if a fire outbreak starting at v 
burns the most vertices of G given an optimal firefighter response. Can the most vulnerable 
vertices in a graph be determined without knowing the optimal firefighter response? Perhaps 
they could then be preemptively defended. From the viewpoint of a network architect, we 
would like to design graphs that are resistant to such attacks. Similar questions can also be 
asked if there are k initial outbreaks of fire. 

Finally, MacGillivray and Wang [B| observed that the firefighter problem can be viewed 
as a one-player game. Suppose that the fire has a choice, too: the fire can only spread 
to d neighbors each time step. This forms a two-player game. What strategy should the 
firefighters use to minimize the number of burned vertices? 
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for 1 < t < 7, 
for 8 < t < T. 
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